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Abstract 

We present, both exactly and approximately, a complete set of mappings between the vacuum (or 
fundamental) leptonic mixing parameters and the effective ones in matter with non-standard neu- 
trino interaction (NSI) effects included. Within the three-flavor neutrino framework and a constant 
matter density profile, a full set of sum rules is established, which enables us to reconstruct the 
moduli of the effective leptonic mixing matrix elements, in terms of the vacuum mixing parameters 
in order to reproduce the neutrino oscillation probabilities for future long-baseline experiments. 
Very compact, but quite accurate, approximate mappings are obtained based on series expansions 
in the neutrino mass hierarchy parameter n = Am^/Am^, the vacuum leptonic mixing parame- 
ter S13 = sin &i3, and the NSI parameters e a n. A detailed numerical analysis about how the NSIs 
affect the smallest leptonic mixing angle #13 , the deviation of the leptonic mixing angle #23 from its 
maximal mixing value, and the transition probabilities useful for future experiments are performed 
using our analytical results. 
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I. INTRODUCTION 



During the past decade, neutrino oscillation experiments have provided us with ver 



convincing evidence that neutrinos are massive and lepton flavors are mixed 



aaa 



is witn very 



101 ] . It opens an important window for searching new physics beyond the Standard 



Model (SM) of particle physics, and has significant cosmological implications. Within the 
framework of three active neutrinos, neutrino masses are the leading mechanism behind 
neutrino orations QflflQ. However, in fn tUr e long-baseline neutrino oscil.ation 
experiments, besides the standard matter effects Ij], H], the possibility of testing non- 
standard neutrino interactions (NSIs) should be opened up. 

Note that, NSIs enter neutrino oscillations at production, propagation, and detection 
processes. In principle, in the case of dimension-6 operators, the corresponding NSI 
parameters are related to the underlying new physics in the form of e ~ (m^/mx) 2 , 
where m^y is the mass of the W boson and mj denotes the new physics energy 
scale. A rough estimate indicates that if new physics appears at the TeV region, the 
magnitude of NSI parameters should not be larger than a few percent, although the 
present experimental upper bounds are still very loose. Due to the interference effects, 
NSIs modify the standard flavor transitions at leading order in e for some typical pro- 



cesses, especial 
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y at a future neutrino 
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or other facilities with hig h-energy beams 
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361 ] . In these cases, 



NSI corrections become particularly relevant, and the experimentally measured values of 
leptonic mixing angles and CP violating quantities are dramatically different from the vac- 
uum ones. In this sense, the combination of standard neutrino oscillations and NSI effects 
in the analyses of future neutrino experiments is not only meaningful but also necessary. 
In addition to NSIs at propagation processes, NSI effects at neutrino sources and detectors 
play a very imp ortant role, since they may induce significant mimicking effects on mixing 
parameters 37[ or bring in distinguishable zero-distance effects for a near detector 38L l39|. 
Here we will only concentrate on NSIs during the phase of neutrino propagation, and a brief 
discussion on how to consistently include the source and detector effects will be implemented 
at the end of Sec. HT1 

Although much attention has been paid on the issue of NSIs according to different neu- 
trino facilities and projects, the previous analytical investigations are either based on a two- 



2 



flavor neutrino framework [40| or an approximation for the three-flavor neutrino oscillation 
probabilities (indeed producing lengthy formulas). There are still lack of analytical relations, 
which can show us the NSI effects on the leptonic mixing parameters in a transparent way. 
Thus, in this paper, we first develop a full set of sum rules, which relate the vacuum leptonic 
mixing matrix elements and their effective counterparts in matter. By solving these sum 



rules, it is straightforward to establish the leptonic mixing matrix, unitarity triangles [41], 
and CP violating effects in matter (see Sec. IHIj) . We then present series expansions of map- 
pings in the mass hierarchy parameter rj = Am^/Am^, the mixing parameter S13 = sin# 13 , 
and the NSI parameters e a p. The NSI corrections to the vacuum mixing parameters can 
be manifested in a distinct way. We hope that the elegant and compact formulas provided 
in this paper could be very helpful for the phenomenological studies of future long-baseline 
neutrino oscillation experiments. 

This work is organized as follows. In Sec. [Til we present the general formulas and no- 
tations, and show how the neutrino oscillation probabilities can be expressed through the 
language of effective mixing parameters. In Sec. IHIl we introduce the sum rules between 
leptonic mixing matrix elements and the effective counterparts in matter, and then derive 
the mappings exactly. The expressions of effective masses in matter, which are necessary 
for the calculation of neutrino oscillation probabilities, are shown in detail in Appendix lAl 
Next, in Sec. lIV| we derive a full set of series expansions of these mappings. We also compare 
our mapping results with the corresponding expressions existing in the literature but with- 
out NSIs, and find that our results are in agreement with previous analyses in the limit of 
vanishing NSIs. Section |V] is devoted to applications of our analytical mappings. Numerical 
illustrations in order to show the validity and reliability of our approximate results are also 
presented. Finally, a brief summary is given in Sec. IVII 

II. THE LANGUAGE OF EFFECTIVE PARAMETERS 

At energy scales /i <C m w , the NSIs involving neutrinos can be described by the effective 
Lagrangian 

£ nsi = ~% E £ S fc-fLvp) (h,Pf) , (1) 

vz f,p 
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where / is a charged lepton or quark, Gf is the Fermi coupling constant, and P = {L, R} is 
a projection operator. The parameters ev^, which are entries of the Hermitian matrix e^ p , 
give the strengths of the NSIs. The magnitudes of the NSI parameters can be constrained 
from neutrino deep inelastic scattering experiments and from elastic v— e scattering in which 
the NSIs would contribute to the determination of sin 2 9w, i.e., the Weinberg angle. The 



latest constraints on have been discussed in Ref. 



42 



43 



44 



45l | , and the most stringent 



bounds are those on ej^ for a = e, /i, r. 

In order to introduce the effective mixing parameters, we start from neutrino oscillations 
in vacuum. The evolution in time of a neutrino state \v(t)) is given by the Schrodinger-like 
equation 

if t \u(t)} = H\u(t)} , (2) 

where H is the Hamiltonian of the system. For neutrinos traveling in vacuum, the Hamil- 
tonian in the ultra-relativistic limit E ^> rrii reads 

if =^Vdiag(0,A 21 ,A 31 )V^ , (3) 
where Ajj = m 2 — m 2 are the neutrino mass-squared differences and E denotes the neu- 



trino energy. In addition, V is the unitary leptonic mixing matrix 



121 ]. which relates the 



mass eigenstates of the three neutrinos (ui, u 2 , v%) to their corresponding flavor eigenstates 



v ™ Ui ' ( 4 ) 



for a = e, fi, r. For simplicity, the sum of Latin indices run over 1, 2, 3 and the sum of Greek 
indices run over e, fi, r throughout this paper, if not otherwise stated. We can define the 
evolution matrix S(t, to) such that 

| u(t)) = S(Mo)K*o)) , S(t ,t ) = l, (5) 

and S(t,to) satisfies the same Schrodinger-like equation (j2J), as \v(t)). The neutrino oscil- 
lation probabilities can be found as P a p = \Sp a (t, to)\ 2 . Using Eq. Q, the elements of the 
evolution matrix are given by 

V ai V; t e~^ , (6) 



where we have identified L = t — to- Thus, the probability of transition from a neutrino 
flavor a to a neutrino flavor (3 is given by 

P a /3 = \S/3a(t, to) | 2 = 

For the setups of future long-baseline neutrino oscillation experiments, the neutrino beams 
inevitably travel through the Earth's mantle, and the charged-current contributions to the 
matter-induced effective potential have to be considered properly. Disregarding the neutral- 
current contributions, the effective Hamiltonian responsible for neutrino propagation in mat- 
ter is given by 

H a j3 = H a p + a (5 ae 5p e + E a p) , (8) 

where the matter parameter a = \f2GpN e arises from coherent forward scattering. Here N e 
denotes the electron number density along the neutrino trajectory in the Earth and the NSI 
parameters e a p are defined as 

f,p e 

with Nf being the number density of a fermion of type /. 

Similar to the vacuum Hamiltonian in Eq. (jHJ), the effective Hamiltonian in matter can 
also be diagonalized through a unitary transformation 

H = ^diag (fhl rh 2 2 , rh 2 3 ) V* , (10) 

where mf denote the effective mass-squared eigenvalues of neutrinos and V is the unitary 
mixing matrix in matter. Note that, in writing down Eq. (110p . we have already taken into 
account the Hermitian property of H. 

In order to write out explicitly the transition probabilities, we assume a constant matter 
density profile, which is actually close to reality for most of the proposed long-baseline 
experiments. Following analogous procedures as shown in Eqs. (H])-([7j), one can then obtain 
the transition probabilities with matter effects included as 

Pa/3 = I S/3 a (t, to) \ 2 = 



E 



(7) 



E 



(11) 
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Comparing Eq. (JTj) with Eq. (jTTj) . we arrive at the conclusion that there is no difference 
between the form of the neutrino oscillation probabilities in vacuum and in matter if we 
replace the vacuum parameters V and mf by the effective parameters V and fhf. The 
mappings between vacuum parameters and the effective ones are sufficient in order to study 
the new physics effects entering future long-baseline neutrino oscillation experiments. The 
key point turns out to be the diagonalization of the effective Hamiltonian H and figuring 
out the explicit relations of the effective parameters. 

In most of the viable models for NSIs, the source and detector effects are simultaneously 
taken into account, despite their magnitude. The language of effective mixing parameters 
can then easily be extended to the case including NSIs at neutrino sources and detectors. 



Now, the NSI parameters at sources and detectors can be defined as 

K) 



17 



31 



46] 



/3=e,fj,,T 



(4\ 



(12) 
(13) 



a=e,/j,,T 

where the superscripts V and 'cf denote source and detector, respectively. The transition 
probabilities are then modified as 1 

2 



a/3 



l + e d y -S(t,t )-(l+e s y 



8a 



^(l+^)^(l+e s ) a5 ^e- i! ^ 



7,5,i 



(14) 



Note that Eq. (114"]) is also suitable to describe neutrino oscillations with a non-unitary mixing 
matrix, i.e., the minimal unitarity violation model [47]]. In the following sections, we will only 
concentrate on NSI effects during propagation processes and establish parameter mappings 
both exactly and approximately. 



III. SUM RULES AND PARAMETER MAPPINGS 

In order to establish analytical relations between the matrix elements of V and those of 
V, we develop a set of sum rules, which enables us to express the products K,iVg* by using 

1 Here we have neglected the normalization factors, which are needed in order to normalize the quantum 
states. 
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V, mf and fhf. Such an approach has partially been employed in Refs. |48|, |49|, [50J in the 
case of three or four-neutrino mixing. Here we will work out the most general form with 
both the standard matter effects and the NSI effects included. 

The first sum rule is just the unitarity conditions, which hold for both V and V, 



E ^ = E = <w • 



(15) 



Inserting Eqs. ([3]) and ( jlOl) into Eq. ([8]), and comparing both sides of this result, it is 
straightforward to obtain the second sum rule 



E rfVjfa = E AnVaiVpi + A a 



(16) 



where we have defined A a p = A (5 ae 5p e + e a p) with A = 2_Ea for simplicity. In order to 
derive a linearly independent sum rule besides the first two, we square both sides of Eq. (Tl6|) 
and obtain the squared relation 

E rtVjfa = E A&VctVZ + E Ary^y + E A ^ (Ary^ + . (17) 

Equations (JT6l) -l fT7l) together with the unitarity condition Eq. ( fTBT) construct a full set of 
linear equations of V^VX (for z = 1,2,3). By solving this set of equations, one will arrive 
at the explicit expressions of VaiV^ straightforwardly. 

In order to be concrete, we reexpress those equations in the following form 



O 



fv n ivA / n \ 



O 



Va2Vp 2 



+ 





A a /3 



A 



A 



v a2 v; 2 



(18) 



where the matrices O, O, and Q are defined by 



/ 



O 



~ 2 ~ 1 ~ 2 

m\ m 2 mg 



~ 4 ~ 4 ~ 4 



1 1 



A 21 A 31 
OA 2 A 2 / 



/ 



Q 





A 21 A 31 



(19) 
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After a lengthy calculation, the solution of Eq. ([TBI can be presented in a very elegant and 
compact form 



V ai V r 



pi 



Aj m Aj n 



A a/3 (m 2 n + m 2 m ) 



(20) 



where i ^ n m, Ay = — mj, and Ay = — rrq — mj. Note that the pairs 
(m, n) = (2, 3), (3, 1), (1, 2) in the right-hand side correspond to % — 1, 2, 3 in the left-hand 
side, respectively. Equation ( 120]) is our main result for the exact analytical mappings. 

The full mappings require the explicit form of m?, which involves the cubic roots of the 
characteristic polynomial of Eq. ([81 . We follow the method given in Ref. 51|, and the 
corresponding roots (or eigenvalues) can be found in Appendix |A] One may worry about 
the m^'s appearing in Eq. (120]) . since it seems that the effective mixing matrix elements rely 
on the absolute effective neutrino masses. However, recalling the solutions in Eqs. (]A4|) - 
(IA6p . it is easy to observe that only the mass-squared differences enter into the expressions 
of rh|, and it guarantees the consistency of our calculations. Obviously, A a p = leads to 



V^V^ = VaiVfo. In the limit e a p — > , Eq . ([20]) red uces to the case of standard matter effects, 



and the main results given in Refs. 



48 



49 



52 



53 



54| will be easily reproduced. These exact 



relations are model independent and do not rely on any specific parametrization, and hence, 
it will be very helpful to systematically study NSIs in future experiments. 

Taking a = (3, the moduli of V ai can be estimated immediately. For the case a ^ (3, 
the sides of leptonic unitarity triangles, which are defined by the orthogonality relations in 
Eq. (|T5]) in the complex plane, are obtained. These unitarity triangles have 18 different sides 
and nine different inner angles, but their areas are all identical to a single rephasing-invariant 
parameter J '/2 defined through j^] 



7,fc 



(21) 



One of the major challenges of future long-baseline neutrino oscillation experiments is to 
measure J7", in order to establish the existence of CP violation in the lepton sector. We can 
also define the counterpart of J in matter as J . Its magnitude is related to the moduli of 
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the effective mixing matrix elements as 

J 2 = l^i 2 !^) 2 !^-! 2 !^! 2 — ^ (i + l^[ 2 |^-I 2 + l^-i 2 !^! 2 

-\V at \ 2 - \% 3 \ 2 - \V ai \ 2 - \V^y . (22) 

As an application, we show, in Appendix [B], the zeroth-order series expansions of | V^>3 1 2 , 
| V^>2 1 2 , and | 1 2 in small parameters, i.e., rj and V e $. In addition, simplified formulas of 
Eq. (1201) and the effective mixing matrix elements for standard matter effects are presented 
in Appendix iBl 

Now, the neutrino oscillation probabilities can be directly obtained with the help of 
Eq. (TTTj) for a realistic experiment. In order to be explicit, we can express the neutrino 
oscillation probabilities in matter as 

P aa = 1 -4^|^.| 2 sin 2 (^f) , (23) 



i>j 



P aP = -4 ]T Re (VaiVpiVajVp^J sin 2 ( ^ j - 8 J J] sin ( ^ j , (24) 

i>j \ J i>j 



where (a,/3) run over (e,/i), (/i,r), and (r, e). For anti-neutrino propagation in matter, we 
can simply recalculate Eqs. (T231 and (|2"4|) through the replacements A — > —A, V ai —> V* t , 
and e a/ 3 — > e* al3 . In general, note that V a i, Ay and J for neutrinos are not identical to 
Vai, Ay and J for anti- neutrinos. At first glance, one may wonder if the information on 
the phases of e a p have been lost, since there is only one parameter J governing the CP- 
violating effects. We stress that, in neglecting the NSIs at sources and detectors, flavor and 
mass eigenstates of neutrinos can always be correlated by using a unitary transformation, 
and hence, we can use one effective rephasing invariant to describe the CP-violating effects in 
neutrino oscillations. For instance, if we ignore the source and detector effects in Eq. (33) of 



Ref . [3l| , the remaining CP-odd terms can be combined together with respect to a common 
oscillating factor, which is consistent with our compact formulas (1231) and (T2"4"|) . 

Although our exact analytical results are very elegant, they do not show how new physics 
affects mixing parameters in a transparent way. From a phenomenological point of view, 
analytically approximate mappings are very useful, since they can reveal the underlying 
correlations between effective mixing parameters and NSI effects, and in particular, show 
which of them are mostly relevant for a given process. In the following section, we will 
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perform a detailed analysis of approximate mappings based on series expansions in small 
mixing parameters and NSI corrections. This method is indeed similar to the analysis of 
series expansion formulas for neutrino oscillation probabilities 56]. 



IV. SERIES EXPANSIONS OF PARAMETER MAPPINGS 

In this section, we proceed to present the series expansion formulas of parameter mappings 
in rj, Sis, and the NSI parameters e a p. For convenience, we adopt the standard parametriza- 
tion and thus the vacuum leptonic mixing matrix V can be parametrized by using three 
mixing angles and one CP violating phase as 



v = o 23 v 6 o 13 vto 



C12C13 S12C13 s 13 e~ lS ^ 

~ S 12C 23 ~~ Ci2S23Si3e l5 C 12 c 23 — Si2S23 s 13 el<S S23 c 13 
\ Sl2 s 23 — C\ 2 C 23 Si 3 e l6 — C 12 S23 — s 12 c 23 s 13 el<5 c 23 c 13 J 



(25) 



where V$ = diag(l, 1, e lS ), and Oy is the orthogonal rotation matrix in the (i, j)-plane with 
dj = cos 8ij and = sin^- (for ij = 12, 13, 23). A global analysis of current experimental 
data yields 0.25 < sin 2 # 12 < 0.37, 0.36 < sin 2 6 23 < 0.67, and sin 2 6b < 0.056 at the 3a 



confidence level, but the CP- violating phase 5 is entirely unrestricted [57J . The best-fit values 
of neutrino mass-squared differences are A21 = 7.65 x 10~ 5 eV 2 and IA31I = 2.4 x 10 -3 eV 2 , 
which indicate that the hierarchy parameter we defined in Sec. [I] is given by 77 = A 21 /A 31 ~ 
±0.032. Present experimental bounds on the NSI parameters e Q/ g show that e MQ (or e afl ) are 
strongly constrained to \e e J < 3.8 x 10~ 4 and —0.05 < £„„ < 0.08 at 90 % confidence 



evel 



30t ] . This is the reason why some of the previous works neglect contributions of e^ a [40| . In 
the following, we will only focus on e er , e^, and e TT contributions, respectively. 

Using a similar notation, one may also define the effective mixing angles #12, #13, #23, an d 
CP violating phase 5 in matter. Then, we can parameterize V in analogy to Eq. (T251) . It is 
straightforward to extract the sines of the mixing angles from Eq. (j25j) using 

Sl3 = |^e3|, Sia= \V e2 \/y/l- \V e3 \ 2 , S 23 =\V fl3 \/\/l-\V e3 \ 2 . (26) 

The effective mixing angles % are obtainable once the moduli of V a p are computed using 
Eq. ( l20l) . and it is not difficult to check that in the limit of vanishing matter effects the 
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effective mixing angles are equal to the vacuum ones. Analytical relations between 9 and 9 
could be very useful and they rely on the perturbation theory that we have employed. 
For our purposes, we first factor out the rotation matrix 23 



H 



A 



31 



23 V 5 ■ M ■ ^/0 2 T 3 



2E 
Agi 
IE 



23 V 5 



V r -diag(A 1 ,A 2 ,A 3 )-V rt 



(27) 



where M is given by 



M = 13 12 • diag (0, n, 1) • Oi 2 Ol 3 + diag [A, 0, Oj + y/o£ ■ e ■ 23 V 5 , 
and A = A/A 3 i. In deriving Eq. (1271) . the following commutative properties are used 



(28) 



vlo 



5^12 



o l2 v ] s , 



V} ■ diag (A, 0, O) = diag (A, 0, o) ■ V} . 



(29) 
(30) 



The diagonalization of M is performed by using perturbation theory, i.e., we write M = 
M (°) + M^ 1 ) + ■ • • , where contains all terms of first order in rj, S\ 3 , and e a p. One finds 



M (0) = diag (A, 0, l) = diag 



(o) x (o) ,(oy 



(31) 



and 



^ijs^ + Aiee r)Si 2 c 12 + Ai efl s 13 e 15 + Ae e 



7]cl 2 + Ae 



V 



ie 7 



/ 



(32) 



with £ ~' denoting the conjugate elements and e Q/ g = (V^ ■ e ■ 23 V$) a p. Since is 
diagonal at zeroth order, we have = I. Then, the first order corrections are given by 



A 



(i) 



M (1) 



(33) 



and 



V; 



(1) 



,(0) ,(0)^ 



(o) 



(34) 



Thus, the effective masses and mixing matrix are given by mf ~ A 31 (A l -°' ) + \\ ) and V ~ 
23 V${V^ + V^), respectively. Finally, inserting Eq. (|32|) into Eqs. (|33|) and ([31]), we arrive 
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at mappings for the effective mass squares 
m\ ~ A 3 i (A + 7]s 2 12 + Ae ee 



m\ ~ A31 



m\ ~ Ac 



rjc\ 2 - Asj 3 ( £ W ~ £ rr) ~ As 23 C 23 (^r + + Ae^ 



1 + Ae, 



L 3 — <-*31 

the effective mixing matrix elements, 
K9 ~ h c 23 e eiM - s 23 e 



) + As 23 c 23 (e MT + e* r ) 



e2 

V, 



A 



er 5 



■si3e lS A(s 23 e efl + c 23 e 
1-A ' 1-A 

„2 



(IT; 



^2 ^ C 23 + S23C23-4 (5 rr - e w ) + S 23 i (s 23 ^T - ci 3 £* T ) 



v; 



/<3 



s 23 + A [s 23 (e w - e rr ) + c 23 £ Mr - S2 3 c 23 (e M r + £^ T ) + S23 {e TT - e w )] 



(35) 
(36) 
(37) 

(38) 

(39) 
(40) 
(41) 



and the effective Jarlskog parameter 



•S13S23 (^f] s sCi 2 c 23 si 2 + Ass+ c f >eti cl 3 \e eiM \ — As$ +< p eT c 23 s 23 \e eT \ \ 
J = (A-l)A 

s 23 ( s 4> eti -4>eT C 2 3 \ e e[j,\ \ter\A 2 — f]S ( j, eti C 12 C 23 Si 2 S 23 \ e efl \ A — T]S^ eT C\ 2 Si 2 c\ 3 \e eT \ A 



(42) 
(A-l)A 

where the a /?'s are the phases associated with the complex NSI parameters e aj 3 and the 
s^Js are the corresponding sine functions. The above mappings can be transferred into 
mappings between mixing angles straightforwardly, i.e., we have approximately §13 ~ \V e3 \, 
§12 — \V e2 \, and §23 — \V^ 3 \. Equations ( l35l) -( )42l) are our main results for the approximate 
analytical mappings. 

Some discussions are in order: 

• In the limit e a p — > 0, it is interesting to observe that our results coincide with the 
mapping results in Ref. 58J] when evaluated at the same order of perturbation theory. 



In fact, expanding Eqs. (27a)- (27c) of Ref. [58| in 77 and S13 (which means that the 
C parameter appearing there equals 1 — A), and retaining terms up to first order in 
these small parameters, we can easily check that they reduce to 

Sl3 



Sl3 
^23 



1 - A 

s 2 3 ? 

V 



Sl2 — ~^C\2 S 12 

A 



(43) 
(44) 
(45) 
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which are in agreement with our Eqs. (l38l) - (|4T!) evaluated at e af 3 = 0. In addition, 
from Eq. (T43l) . it is worth noticing that the mixing parameter S13 is strongly modified 
by matter effects when < A 7^ 1, otherwise the resonance A = 1 is at work and the 
mapping procedure adopted in this paper is not valid. 

As shown in Eq. (129ft . the orthogonal matrix 23 commutes with the standard matter 
potential. Hence, in the case of vanishing 9 13 and NSIs, Eq. (jHJ) can be rewritten as 



A31 

H = — 93 
2E 



12 • diag (0, rj, 1) • Of 2 + diag (A, 0, o)] 2 T 3 , (46) 



where the CP violating phase 5 loses its meaning and does not appear. An evident 
conclusion deduced from Eq. (146j) is that, if #13 = 0, the standard matter effects 
only contribute to the mixing angle #i 2 . Both $13 and # 2 3 as well as A31 will keep 
their vacuum values in matter. However, when the NSIs are taken into account, the 
situation will be quite different. A non-zero #13 will emerge in general, # 2 3 will deviate 
from its maximal value 7r/4, and a non-trivial CP violating phase 5 may also exist. 

As already noticed above, the mapping for 9 12 shows an unphysical divergence for A — > 
or A — > 1, and the vacuum results cannot be reproduced, a well-known consequence 
of the perturbative approach adopted in the mapping procedure. Thus, degenerate 
perturbation theory should be elaborated on around these two singularities. 

Except from e efl and e er , it can also be very clearly seen that contributions to §13 
from all the other NSI parameters are all suppressed. Since the present experimental 



bound on e eii is rather stringent [43J, we conclude that e er is the most significant NSI 
parameter to be taken into account for # 13 . As for # 23 , NSI corrections are relatively 
mild unless very high-energy regions are considered. 

The matrix elements |V^3|, | 1 , |V^i|, and |Vrt| are not modified by e eT . 



V. APPLICATIONS 

We now proceed to numerically illustrate (using normal mass hierarchy, i.e., rj ~ +0.032) 
the NSI corrections to the leptonic mixing parameters and the neutrino oscillations probabil- 
ities based on our model independent results obtained in Sees. IIHI and IIVI We first consider 
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the exact results for mixing matrix elements, and combinations of them, as obtained ap- 
plying Eq. (I20p . This is an important point, since both mixing angles and the Jarlskog 
parameter depend on the modified behavior of the effective matrix elements [see Eqs. (12 2p 
and 0261) ] due to matter effects and non-standard physics. We then focus on the importance 
of the NSI effects o n 9i 3 and #23 in correcting their tri-bimaximal mixing values (i.e., #13 = 



and 823 = 7r/4) 5jJ|60||. Finally, in order to show the goodness of our approximate results for 



the effective mixing angles, we compare the exact v e — > z/ M , v e — > i/ T , and — > v T oscillation 
probabilities obtained including the NSI effects with those derived using Eqs. (l38l) - (f4"T!) for 
the effective mixing angles. 



A. NSI corrections to the leptonic mixing matrix 

The effective leptonic mixing matrix can be reconstructed directly from Eq. (1201) . For 
a = (3, we obtain the expressions of the matrix elements \V a i\, whereas for a 7^ /3, we obtain 
the sides of unitarity triangles V Q jVg* . The numerical results are presented in Fig. [TJ In our 
numerical calculations, we take the central values of the neutrino mass-squared differences 
(A 2 i = 7.65 x 10~ 5 eV 2 , A 3i = 2.4 x 10~ 3 eV 2 ) and the leptonic mixing angles (9 12 = 33.5°, 
623 = 45°, #13 = 0) obtained in a global analysis of the presently available neutrino oscillation 
as an example, we choose e eT as the only non-vanishing NSI parameter. 
For comparison, we also show the results without including NSIs. In particular, for higher 
neutrino beam energies, the NSI corrections are remarkable. 

First, we can observe that the energy dependence of the matrix elements can be easily read 
off from the approximate relations in Eqs. (!38l) - (f4TI) . since the matter parameter A ~ E? 
Thus, for example, the fact that \ V e 2\ is predicted to decrease with increasing neutrino energy 
is confirmed in the first panel of the fourth row in Fig. [TJ Moreover, the singularity for matrix 
elements around E ~ 10 GeV (see panels in the first and third rows) clearly corresponds to 
the resonance at A ~ 1, which can be understood investigating the perturbative results in 
Eqs. (l38|) - (f4TI) . Note that, for anti-neutrinos, since the sign in front of A is negative, such a 
singularity should not appear. In addition, the relation |KiV^ | — iKi^l holds quite well, 
which is an obvious consequence of the fi — r symmetry in the genuine neutrino mass matrix. 

2 The matrix elements that are not quoted in Eqs. ([38|) - (f4T|) can be obtained using unitarity relations. 



14 



10 



10" 



IS- 



10" 





NSI: \V^V* m \ 





NSI: |V ei V;*| 



r 

si: |y Ml | 2 ! 



I 

\ I 
N i 



10" 



10" 



NSI: \V el \ 2 


T= 




\ 

\ 




\ 




\ 











NSI: \V^\ 2 



I 

\ / 
\ / 



r 



\ I : 
\ / 

NSI: 



SI: \V T , 



NSI: \V Tl \ 2 



I 

\ / 
N / 



10° 10 1 



10 2 10° 



10' 

E (GeV) 



10 2 10° 



10 



1(T 



FIG. 1: Illustrative plots for the effective matrix elements (first and second rows) and \V a i\ 2 

(third and fourth rows) as a function of the neutrino energy E. Here the solid, dashed, and dotted 
curves correspond to i = 1,2,3, respectively. The first and third rows show the results without 
including NSIs (labeled SI), while the second and fourth rows are those including NSIs. We use the 
representative value Re(e er ) = Im(e er ) = 0.02, with all other E a p being zero. 
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This can be seen comparing the first and second panels of the second row for any value of 
the index i. In the same panels, at energies E ~ 4 GeV, V e iVZ and V e iV*i (i = 1, 2, 3) 
are comparable to each other, and thus, the unitarity triangle built with these sides takes 
a nearly equilateral form with three nearly degenerate inner angles. Such an equilateral 
form is destroyed when increasing the energy. Similarly, IV^iV^-J is rather stable against 
matter corrections and NSI effects, which also reflects the stabilization of 823 . As for the 
matrix elements, |V^i| and \ V T \ \ are not sensitive to e eT , which is also in agreement with our 
approximate mappings. 



B. NSI corrections to the mixing angles 

A crucial goal of future neutrino facilities is to measure the smallest leptonic mixing 
angle # 13 in order to extract information on leptonic CP violation. However, it has been 
pointed out that NSIs may play a very important role for mimicking effects on #13 and 



leptonic CP violation, especially in the case of a small 6*13 6JJ. It is then quite important 
to analyze in detail these effects in order to be able to disentangle genuine #13 effects from 
new physics-induced ones at future neutrino facilities. 

On the other hand, the question of whether the leptonic mixing angle 823 is exactly 
maximal or not is quite relevant, especially from the model builders' point of view: in 
fact, many models presented in the literature predict 823 being (almost) maximal and the 
understanding of the flavor problem strongly relies on the knowledge on the value of 823 
to be as accurate as possible. Thus, it is very important to investigate the possible NSI 
corrections to #13 and the maximal mixing pattern in the fx — t sector. 

According to Eqs. (|38p - (l4ip . the most relevant NSI parameter for # 13 is e eT (since the 
upper bound on e etl is rather stringent), whereas for 6*23 e^ T , and e TT contribute. Notice 
that, in the latter case for maximal mixing (#23 = 45°), a typical feature is that the vacuum 
Hamiltonian takes on a fi — r symmetric form, namely, H is invariant under the exchange 
of fx and r indices. Hence, if NSIs possess a similar fx — t symmetric form (i.e., e efJi = e eT 
and £ MjU = e TT ) , the fx — t symmetry exists in the effective Hamiltonian H, and the effective 
mixing angle 823 will not be affected by matter effects. As a consequence, e^ T itself does not 
contribute to 823 if all the other NSI parameters are zero. 

In the upper plots of Fig. [2j we show the non- vanishing 6*13 generated by the NSIs [com- 
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FIG. 2: Dependence of 'the mixing angles #13 and 623 on the NSI parameters for three representative 
values of the neutrino energy: E = 1,5, and 30 (35) GeV, which roughly correspond to the v e and 
Vn mean energies at a 50 GeV neutrino factory. The vacuum value of O13 is fixed to be zero, whereas 
we assume maximal mixing for 623 ■ In each plot, the darker the region, the larger the deviation of 
9\3 and 623 from their vacuum values. Only the labeled NSI parameters are non-vanishing in each 
plot. 



puted using our exact formula given in Eq. ([20]) ]. One can observe that 613 is quite sensitive 
to e er . In the case of E = 30 GeV, #13 may acquire a very sizable value close to 90°. This 
is due to the reordering of the eigenvalues rhi and m 3 when E > 10 GeV. If we keep the 
order of eigenvalues in the form of diag(m^, m|, m|), a shift of it/2 has to be added to #13. 
In the lower plots of Fig. [2, we show the corrections from £ w and e TT to 623, assuming the 
vacuum value 823 = ir/A. Our numerical results indicate that there are no sizable corrections 
to the mixing angle 823, and even in the high-energy region, 823 should not deviate from its 
maximal value by more than a few degrees. We also find that the e MM — e TT contributions 
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FIG. 3: Neutrino energy dependence of the effective mixing angles (19 13 and #23 fib/id curves 
correspond to exact numerical results, whereas the dashed ones are computed using our approximate 
mappings. The non-vanishing NSI parameters have been labeled in each plot. 



are symmetric with respect to the e m = e rr axis up to a minus sign, since an addition to 
one of the parameters could as well be made to the other one. 

It is interesting to observe that the main features of the previous exact results can also 
be captured from our approximate mappings in Eqs. (j38j) - (HIl) . To illustrate this point, we 
show the dependence of #13 and #23 on the neutrino energy in Fig. [3J for different values of 
the relevant NSI parameters, according to our foregoing discussions. Solid curves correspond 
to exact results obtained using Eq. (1201 . whereas dashed ones represent our perturbative 
mappings. In the first row, we can appreciate how the dependence of #13 on e eT is well 
described by our perturbative result in Eq. ( 1381) . unless e eT assumes a very large value, close 



to its upper bound 



301 ]. In addition, notice that the increase of 6 I 13 corresponds to reordering 



the eigenvalues for energies around 10 GeV. In the second row, we analyze the behavior of 
#23, for different values of the relevant parameters and e TT . In particular, in the first and 
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second panels, we choose only one of them being different from zero (and equal to 0.01), 
whereas in the last one, we allow both of them to assume larger values (e w = 0.08 and 
e TT = 0.4). The agreement between our calculation and the exact evaluation of 6*23 is quite 
good, also in predicting the location of the resonance. 

Finally, we comment on the fact that 9 12 is dramatically suppressed by matter effects, 
as shown in Eq. (|38|) . However, since long-baseline neutrino oscillation experiments are not 
very sensitive to this angle, we will not perform a detailed analysis here. The conclusions 
made above about the dependence of the effective angles on the NSI parameters apply as well 
to the case of anti-neutrinos in matter and we will not perform such a redundant analysis 
here. 



C. NSI corrections to the neutrino oscillation probabilities 

Until now we have described the relevant features of new physics effects on the matrix 
elements of the leptonic mixing matrix. In what follows, we study the dependence of the 
transition probabilities on the NSI parameters, for different choices of neutrino energies 



and baselines. In particular, we focus on the golden channel v { 



asymmetries derived from it, the silver channel v e — > v T 63|, [64j, and the so-called discovery 



62| and the CP 



channel — > u T , which is thought to be the best channel for searching for new physics [65]. 
We also show how the relevant features of the transition probabilities are well reproduced 
computing them by inserting Eqs. (T38|) - (l4TT) into Eq. (|TT|) . 

In Fig. HI we show the transition probability P{y e — ► v^) as a function of the neutrino 
energy for three different baseline setups: L = 700 km (around the scope of MINOS 



and OPERA 66j), L = 3000 km, and L = 7000 km (for the two detector setup of a 
neutrino factory). The input parameters are the same as those in Fig. [TJ In each panel, 
the solid curves denote the exact numerical results, the dashed curves correspond to results 
derived from our approximate mappings and, to highlight the effects of the NSI parameters, 
the dotted curves represent the probability without including NSIs. We can observe that 
our approximate mappings given in Sec. [IV] agree with the exact numerical results to an 
extremely good precision. Similar to the plots of the mixing parameters, a singularity exists 
around E ~ 10 GeV due to the limitation of non-degenerate perturbation theory that we 
have elaborated. For smaller #13 , the probability is more sensitive to the NSI effects, and 
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FIG. 4: Neutrino oscillation probabilities for the u e — > channel as a junction of the neutrino 
energy E. The baseline lengths and values of S13 have been labeled in each plot. Here, we set 
5 = tt/2 and only e eT = 0.01 is allowed to be non-vanishing. The solid curves denote the exact 
numerical results. The dashed curves correspond to results derived from our approximate mappings, 
and for comparison, the dotted curves are shown to illustrate probabilities without including NSIs. 



thus, longer baseline lengths are more favored for the purpose of searching for new physics 
effects. 

The experimentally measured CP asymmetry in the golden channel, which is usually 
defined as 

a P(z/ e -> iy) - P(V e -> gg) (A7 , 

is illustrated in Fig. [5] for the same baseline setup. Again, our approximate mappings 
are valid in a large range of beam energies. At higher energies, the CP asymmetries are 
dramatically affected by NSIs, i.e., the values of Acp, which are calculated without taken 
into account NSI effects, may go to divergent directions. 

In Figs. [6] and we repeat the same exercise on the neutrino oscillation probabilities 
and CP asymmetries, but instead as a function of the baseline length and for two fixed 
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FIG. 5: CP asymmetry Ac p derived from the v e — > channel. The values of the mixing parameters 
as well as those of the baselines and neutrino energies are the same as in Fig. The solid 
curves denote the exact numerical results, the dashed curves correspond to results derived from our 
approximate mappings, and the dotted curves show the probabilities without including NSIs. 



value of the neutrino energy E = 5 GeV and E = 30 GeV. It can be clearly seen that 
for smaller 8 13 and lower beam energy, new physics effects play a significant role around 
L ~ 3000 km, which sheds some light on future beta beam experiments. For higher energy 
experiments, i.e., a neutrino factory, a far detector with relatively longer baseline length 
should be important to constrain NSI effects. In both figures, one can appreciate how the 
probabilities computed using our approximations for the effective mixing angles are in very 
good agreement with the exact results. 

Finally, we illustrate the application of our analytical expressions for the v e — > v T and 
— > v T channels in Fig. [SJ For comparison, we also show the maximal NSI corrections by 
setting all the NSI parameters at their upper bounds given in Ref. 



30J. One can observe 

that NSI corrections to these two channels are not remarkable if the corresponding e a ps are 
chosen to be a few percent. However, increasing the NSI parameters, the NSI effects become 
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FIG. 6: Transition probability P(y e — > v^) as a function of the baseline length L. Here, e er = 0.01 
and 5 = it/2 are adopted, and the neutrino beam energies have been labeled in each plot. The 
solid and dotted curves denote the exact numerical results with and without NSIs, respectively. 
Probabilities calculated using our approximate mappings are shown as dashed curves. 

more significant, in particular for the v e — > v T channel. The upper plots in Fig. [8] indicate 
that our approximate mappings are not quite valid for relatively longer baseline lengths. This 
is due to the fact that our expansions are performed according to small £ a pS and cannot be 
extended to the regions of sizable NSI parameters. As discussed in the introduction, if NSIs 
are related to some underlying new physics, they should be attributed to next-to-leading 
order effects and not deviate much from zero. In this sense, our approximate mappings 
are quite realistic and should be very helpful for both phenomenological studies and model 
buildings. 

Since the analyses above certainly depend on the input NSI parameters, they mainly 
serve as illustrations. However, our analytical results are model independent. Thus, they 
are hoped to be very useful for a general study of NSI effects in future experiments. The 
transparent mappings also manifest the underlying correlations between leptonic mixing 
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FIG. 7: CP asymmetries as a function of the baseline lengths. Here, e er = 0.01 and the neutrino 
beam energies have been shown in each plot. The solid and dotted curves denote the exact numer- 
ical results with and without NSIs, respectively. Dashed curves denote Acp calculated using the 
approximate mappings. 



parameters and NSIs in a very legible way. 



VI. SUMMARY 



In this work, we have developed both exact and approximate mappings between the 
leptonic mixing matrix in vacuum and in matter in the presence of NSIs. A full set of 
sum rules between fundamental mixing parameters and the corresponding effective ones in 
matter have been derived. By using these sum rules, exact and model independent analytical 
mappings between the mixing matrix elements V a i and V a i have been established, and in 
turn using these mappings, the moduli of the mixing matrix elements and the sides of 
unitarity triangles can be immediately figured out. Besides the exact expressions for the 
mixing parameters, we have also derived approximate parameter mappings based on series 
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FIG. 8: Transition probabilities for the v e 



u T and v.. 



v T channels as a function of the 



baseline length L. Different setups of N SI parameters are considered: (a) for the u e — ► v T channel, 



we choose e e 



0.01, and for the za, 



v T channel, we choose and £ 



/'/' 



0.01. 



Here, thick solid and dashed curves correspond to the exact numerical results and our approximate 



M] with 



mappings, respectively; (b ) NSI parameters are at their upper bounds computed in Ref. 
thin solid and dashed curves corresponding to the exact numerical results and our approximate 
mappings, respectively. Dotted curves denote the numerical results without including NSIs, and 
they are unique in each channel. 



expansions in the small parameters 77, S13, and e a p- We have then performed a detailed 
numerical analysis of the application and validity of our parameter mappings. In particular, 
we have concentrated on the mimicking effects of NSIs on the mixing angle #13 and on the 
deviation of the mixing angle 623 from maximal mixing. Furthermore, we have studied in 
detail how the e a ps affect the transition probabilities of the u e — > z/^, v e — > u T , and — > v T 
channels. We have found that the exact parameter mappings are very useful in obtaining 
exact results for the mixing parameters and transition probabilities, and our perturbative 
parameter mappings also describe quite well all the relevant features of these quantities. 
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Note that our analytical analysis is independent of any specific model or assumptions on the 
configuration of NSI parameters. In conclusion, the outstanding feature of our parameter 
mappings is that they reveal the underlying correlations between NSI effects and neutrino 
mixing parameters in a highly straightforward way, and they are very practical and useful 
for the study of NSIs in future long-baseline neutrino oscillation experiments. It also makes 
sense to note that the calculation procedures we have employed in the present work can also 



be applied to the picture of non-unitary leptonic mixing 
elsewhere. 



471 ]. which will be elaborated on 
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APPENDIX A: CALCULATIONS OF EFFECTIVE MASSES 



To calculate the explicit expressions of rhi, the cubic roots of the characteristic polynomial 



of Eq. (jSJ) are involved. We follow the method given in Ref. 
elementary symmetric polynomials 



67 



68 



51] and define the so-called 
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It is easy to check that the relations c 2 = — J2i ™l /(^E), c\ = Yli<j mj-fh^/^E) 2 , and 
Co = — Yli ™i /(2E) 3 are satisfied. By incorporating the definitions above, the mass squared 
eigenvalues can be computed as 
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where p = c\ — 3ci and q = — 27c /2 — c 3 , + 9cxc 2 /2. As a natural consequence, the effective 
mass eigenvalues in matter are only related with the neutrino mass squared differences but 
not the absolute neutrino masses. 

As an example, we consider the case of vanishing NSI. From Eqs. f)Alj) - (lA3j) . one can 
directly write down 
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Substituting Eqs. (!A7j) - (lA9p into Eqs. (1A4h - (1A6[) . the matter corrected eigenvalues given in 



Refs. 
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70| | can be reproduced straightforwardly. One may also check that A = leads to 



the limit rhi = rrii 
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APPENDIX B: FORMULAS FOR EFFECTIVE MIXING MATRIX ELEMENTS 



In the limit of small parameters, i.e., 77 — and V e3 — > 0, one can use our main result for 
the exact analytical parameter mappings Eq. f[2"Ul) to derive zeroth-order series expansion 
formulas for the modulus squares of the mixing matrix elements V e3 , V e2 , and V^ 3 . The 
results are given by 
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-Ae m (m 2 + m 2 2 ) + 2AA 31 [e^\V^\ 2 + Re(e e ^V e3 V* 3 ) + Re(e /1T y T3 V^)] } 



(B3) 



which are valid to all orders in the NSI parameters. In addition, for standard matter effects, 
i.e., without NSI effects, and for any r] and T4 3 , we can derive the corresponding formula to 
Eq. ([20]). The result is 
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In the specific cases of V e3 , V e 2, and V^, we obtain 
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which are valid to all orders in the small parameters 77 and V e3 . 
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